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ABSTRACT 



Aims. To study scalar tensor theories of gravity with power law scalar field potentials as cosmological models for accelerating universe, 
using cosmic clocks. 

Methods. Scalar-tensor quintessence models can be constrained by identifying suitable cosmic clocks which allow to select confidence 
regions for cosmological parameters. In particular, we constrain the characterizing parameters of non-minimally coupled scalar- tensor 
cosmological models which admit exact solutions of the Einstein field equations. Lookback time to galaxy clusters at low intermediate, 
and high redshifts is considered. The high redshift time-scale problem is also discussed in order to select other cosmic clocks such as 
quasars. 

Results. The presented models seem to work in all the regimes considered: the main feature of this approach is the fact that cosmic 
clocks are completely independent of each other, so that, in principle, it is possible to avoid biases due to primary, secondary and so 
on indicators in the cosmic distance ladder. In fact, we have used different methods to test the models at low, intermediate and high 
redshift by different indicators: this seems to confirm independently the proposed dark energy models. 

Key words, cosmology: theory - cosmology: quintessence -Noether symmetries-Scalar tensor theories 



1. Introduction 

An increasing harvest of observational data seems indicate that 
^ 70% of the present day energy density of the universe 
is dominated by a mysterious "dark energy" component, de- 
scribed in the simplest way using the well known cosmolog- 
ical constant A ( Perlmutter et al., 1997 Perlm utter et al., 19991 
[Riess et al., 1998, ,Riess, 200 ) and explains the accelerated ex- 
pansion of the observed Universe, firstly deduced by luminosity 
distance measurements. However, even though the presence of a 
dark energy component is appealing in order to fit observational 
results with theoretical predictions, its fundamental nature still 
remains a completely open question. 

Although several models describing the dark energy compo- 
nent have been proposed in the past few years, one of the first 
physical realizations of quintessence was a cosmological scalar 
field, which dynamically induces a repulsive gravitational force, 
causing an accelerated expansion of the Universe. 

The existence of such a large proportion of dark energy in 
the universe presents a large number of theoretical problems. 
Firstly, why do we observe the universe at exactly the time in 
its history when the vacuum energy dominates over the matter 
(this is known as the cosmic coincidence problem). The second 
issue, which can be thought of as a fine tuning problem, arises 
from the fact that if the vacuum energy is constant, like in the 
pure cosmological constant scenario, then at the beginning of 
the radiation era the energy density of the scalar field should 
have been vanishingly small with respect to the radiation and 
matter component. This poses the problem, that in order to ex- 
plain the inflationary behaviour of the early universe and the late 
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time dark energy dominated regime, the vacuum energy should 
evolve and cannot simply be constant. 

A recent work has demonstrated that the fine-tuning problem 
can be alleviated by selecting a subclass of quintessence models, 
which admit a tracking behaviour ( |Steinhardt et al, 1999| ), and 
in fact, to a large extent, the study of scalar field quintessence 
cosmology is often limited to such a subset of solutions. In 
scalar field quintessence, the existence condition for a tracker 
solution provides a sort of selection rule for the potential V{(p) 
(see ( |Rubano et al., 20()4| l for a critical treatment of this ques- 
tion), which should somehow arise from a high energy physics 
mechanism (the so called model building problem). Also adopt- 
ing a phenomenological point of view, where the functional 
form of the potential V{4>) can be determined from observational 
cosmological functions, for example the luminosity distance, 
we still cannot avoid a number of problems. For example, an 
attempt to reconstruct the potential from observational data 
(and also fitting the existing data with a linear equation of state) 
shows that a violation of the weak energy condition (WEC) is 
not completely excluded ( [Caldwell et al., 2003), and this would 
imply a superquintessence regime, during which < — 1 
(phantom regime). However it turns out that, assuming a dark 
energy component with an arbitrary scalar field Lagrangian, 
the transition from regimes with > -1 to those with 
W0 < -1 (i.e. crossing the so called phantom divide) could 
be physically impossible since they are either described by a 
discrete set of trajectories in the phase space or are unstable 
( |Hu, 2005] [Vikman, 2005[ ). These shortcomings have been 
recently overcome by considering the unified phantom cosmol- 
ogy (Nojiri & Odintsov, 2005 ) which, by taking into account a 
generalised scalar field kinetic sector, allows to achieve models 
with natural transitions between inflation, dark matter, and dark 
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energy regimes. Moreover, in recent works, a dark energy com- 
ponent has been modeled also in the framework of scalar tensor 
theories of gravity, also called extended quintessence (see for 



Dvali et al., 2000) Chiba, 1999| 
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Uzan, 1998 


Capozziello, 2002 Capozziello et al., 2003 


Nojiri & Oc 


intsov, 2003 |Vollick, 2003| |Meng & Wang, 2003 


Allemandi et al., 2004|i. Actually it turns out that they are 



compatible with a peculiar equation of state w < -1, and 
provide a possible link to the issues of non-Newtonian gravity 
( |Fujii, 2000) l. In such theoretical background the accelerated 
expansion of the universe results in an observational effect of 
a non-standard gravitational action. In extended quintessence 
cosmologies (EQ) the scalar field is coupled to the Ricci scalar, 
R, in the Lagragian density of the theory: the standard term 
167rG, R is replaced by \6nF{(fi) R, where F((p) is a function 
of the scalar field, and G, is the bare gravitational constant, 
generally different to the Newtonian constant Gn measured in 
Cavendish-type experiments ( |Boisseau et al., 2000| ). Of course, 
the coupling is not arbitrary, but it is subjected to several con- 
straints, mainly arising from the time variation of the constants 
of nature ( |Riazuelo & Uzan, 2002| l. In EQ models, a scalar field 
has indeed a double role: it determines at any time the effective 
gravitational constant and contributes to the dark energy density, 
allowing some different features with respect to the minimally 
coupled case (|Riazuelo & Uzan, 2002]l. Actually, while in the 
framework of the minimally coupled theory we have to deal with 
a fully relativistic component, which becomes homogeneous on 
scales smaller than the horizon, so that standard quintessence 
cannot cluster on such scales, in the context of nonminimally 
coupled quintessence theories the situation changes, and the 
scalar field density perturbations behave like the perturbations 
of the dominant component at any time, demonstrating in the so 
called gravitational dragging ( |Perrotta et al., 2000| . 

In this work, we focus our attention on the effect that dark en- 
ergy has on the background evolution (through the t(z) relation) 
in the framework of some scalar tensor quintessence models, for 
which exact solutions of the field equations are known. In par- 
ticular we show that an accurate determination of the age of the 
universe together with new age determinations of cosmic clocks 
can be used to produce new strict constraints on these dark en- 
ergy models, by constructing the time-redshift (f - z) relation, 
and comparing the theoretical predictions with the observational 
data. 

Although discrepancies between age determinations have 
long plagued cosmology, the situation has changed dramatically 
in recent years: type la supernova measurements, the acoustic 
peaks in the CMB anisotropies ( Spergel et al., 2003 1, and so on, 
are all consistent with an age to - 14+ IGyr. Recently Krauss & 



Chaboyer ( Krauss & Chaboyer, 2003 ) have provided constraints 



on the equation of state of the dark energy by using new globular 
cluster age estimates. 

Furthermore, Jimenez & al. proposed a non parametric 
measurements of the time dependence of w{z), based on the 
relative ages of stellar populations pimenez & Loeb, 2002[ 
[Jimenez et al., 2003) . It is therefore rather timely to investigate 
the implications of new age measurements within in the frame- 
work of quintessence cosmology. 

It turns out that such (f - z) relations are strongly varying 
functions of the equation of state w(z), and could also be very 
useful for breaking the degeneracies that arise in other observa- 
tional tests. 

As a first step toward this goal, we study the cosmo- 
logical implications arising from the existence of the quasar 



APM 08279+5255 and three extremely red radio galaxies at 
z = 1.175 (3C65), z = 1.55 (53W091) and z = 1.43 
(53W069) with a minimal stellar age of 4.0 Gyr, 3.0 Gyr 
and 4.0 Gyr, respectively, extending the analysis performed 
in ( [Lima & Alcaniz, 20051 [Alcan itz et al., 20 03) to a scalar 



tensor field quintessential model (Capozziello & de Ritis, 1994 
[Demianski et al., 2006| [Capozziello et al, 1996f ! 

As a further step, we also consider the lookback time 
to distant objects, which is observationally estimated as the 
difference between the present day age of the universe and 
the age of a given object at redshift z, akeady used in 
(Capozziello et al., 2004) to constrain dark energy models. Such 
an estimate is possible if the object is a galaxy observed in more 
than one photometric band, since its color is determined by its 
age as a consequence of stellar evolution. 

It is thus possible to get an estimate of the galaxy age by 
measuring its magnitude in different bands and then using stellar 
evolutionary codes. It is worth noting, however, that the estimate 
of the age of a single galaxy may be affected by systematic errors 
which are difficult to control. 

It turns out that this problem can be overcome by consid- 
ering a sample of galaxies belonging to the same cluster In 
this way, by averaging the age estimates of all the galaxies, 
one obtains an estimate of the cluster age and thereby reduc- 
ing the systematic errors. Such a method was first proposed by 
Dalai & al. (Dalai etal., 2001) and then used by Ferreras & al. 
( [Ferreras et al., 2003 1 to test a class of models where a scalar 
field is coupled to the matter term, giving rise to a particular 
quintessence scheme. We improve here this analysis by using 
a different cluster sample fAndreon et al, 2004| [Andreon, 2003| ) 
and testing a scalar tensor quintessence model. Moreover, we 
add a further constraint to better test the dark energy models and 
assume that the age of the universe for each model is in agree- 
ment with recent estimates. Note that this is not equivalent to the 
lookback time method as we will discuss below. 

The layout of the paper is as follows. In Sect. 2, we briefly 
present the class of cosmological models which we are going to 
consider, defining also the main quantities which we need for the 
lookback time test. Sects. 3 and 4 are devoted to the discussion 
of cosmic clocks at low, intermediate and high redshifts whose 
observational data are used to test the theoretical background 
model. Finally we summarize and draw conclusions in Sect. V. 
The lookback time method is outlined in the Appendix A, refer- 
ring to (CapozzieUo et al., 2004| ) for a detailed exposition. 



2. The model 

The action for a scalar-tensor theory, where a generic 
quintessence scalar field is non-minimally coupled with grav- 
ity, and a minimal coupling between matter and the quintessence 
field is assumed, is: 



v(0) + x« 



(1) 



where F(0), and y(<^) are two generic functions, representing 
the coupling with geometry and the potential respectively, R 

is the curvature scalar, -g'^^(p^f,(p^v is the kinetic energy of the 

quintessence field <p and describes the standard matter con- 
tent. In units 87rG = fi = c - 1 and signature (+, -, -, -), we 
recover the standard gravity for F = -1/2, while the effective 

gravitational coupling is Gg/f - -- — . Chosing a spatially flat 
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Friedman-Robertson-Walker metric in Eq. ([TJ, it is possible to 
obtain the pointlike Lagrangian 



X. - GFad"^ + 6F'(f>a^d + a^p^ - D , 



(2) 



where a is the expansion parameter and D is the initial dust- 
matter density. Here prime denotes derivative with respect to 
and dot the derivative with respect to cosmic time. The dynami- 
cal equations, derived from Eq. are 



" ^'2f[-^-I 



2H + 3H^ = —p^, 



(3) 



(4) 



where the pressure and the energy density of the ^-field are given 
by 



p^ = - V(<f>) - 2(F + 2HF) , 



+ V((p) + 6HF , 



(5) 



(6) 



and p„, is the standard matter-energy density. Considering Eq. 
(O, the variation with respect to (p gives the Klein-Gordon equa- 
tion 



^ + 3H^ + 12WF' + 6HF' + V' ^0. 



2.0.1 . The case of power law potentials 



(7) 



Recently it has been shown that it is possible to determine the 
structure of a scalar tensor theory, without choosing any specific 
theory a priori, but instead reconstructing the scalar field poten- 
tial and the functional form of the scalar-gravity coupling from 
two observable cosmological functions: the luminosity distance 
and the linear density perturbation in the dustlike matter com- 
ponent as functions of redshift. Actually the most part of works 
where scalar-tensor theories were considered as a model for a 
variable A-term followed either a reconstruction point of view, 
either a traditional approach, with some special choices of the 
scalar field potential and the coupling. However it is possible 
also a third road to determine the structure of a scalar tensor 
theory, requesting some general and physically based properties, 
from which it is possible to select the functional form of the cou- 
pling and the potential. An instance of such a procedure has been 
proposed in (Demians ki et al., 1991[ Capozziello et al, 1996 1: it 
turns out that requiring the existence of a Noether symmetry for 
the action in Eq. ([Til, it is possible not only to select several ana- 
lytical forms both for F{(f)) and V{(f)), but also obtain exact solu- 
tions for the dynamical system ( I3I4I7I ). In this section we analyze 
a wide class of theories derived from such a Noether symmetry 
approach, which show power law couplings and potentials, and 
admit a tracker behaviour. Let us summarize the basic results, 
referring to (I Demianski et al., 2006| [Demianski et al., 2007[ ) for 
a detailed exposition of the method. It turns out that the Noether 
symmetry exists only when 



V = Vo(Fi<p)y^'\ 

where Vq is a constant and 
3(s+l) 



(8) 



Pis) 



where s labels the class of such Lagrangians which admit a 
Noether symmetry. A possible choice of F{(p) is 



F = ^(s)(4, + <Pif , 
where 

{2s + 3f 



48(,s + \){s + 2) 



(10) 



(11) 



and (pi is an integration constant. The general solution corre- 
sponding to such a potential and coupling is: 



a(t) = A{s) \B{s)t- + —1 f , 

4,{t) ^ C(s)\-^B{s)t- + —\ r— , 

where A(s), B(s), C(s), yis) andxis) are given by 
A(s) = (x(s)Y 

B(s) 



3r(s) 

(s + 3)^oY^> (^ + 3)' 



3r(^) 

C(s) = (x(s)) 
and 



s + 6 



_<3+2s) 

(2£i3) /(5H-3)So\ 



3r(^) 



y(s) = 



2s + 3 



X(s) = - 



12(s + l)(s + 2) ' 
2s 



2s + 3' 



(12) 
(13) 

(14) 
(15) 
(16) 

(17) 
(18) 



where D is the matter density constant, So is an integration con- 
stant resulting from the Noether symmetry, and Vq is the constant 
which determines the scale of the potential. Even if these con- 
stants are not directly measurable, they can be rewritten in terms 
of cosmological observables like Hq, Q.^ etc, as in the following: 



D 



1 

Ms) 



B(s) 



(19) 



<^Hs + 3)' 



+3 (s + 6) X 



{(Ho-2)s^ + 3(Ho-2)s-3)y(s)i 



^ \ A -1-3 



(^ + m^) 



(20) 



Here we are following the procedure used in 
( IDemianski et al., 2006| ), taking the age of the universe, to, 
as a unit of time. Because of our choice of time unit, the ex- 
pansion rate H(t) is dimensionless, so that our Hubble constant 
is not (numerically) the same as the Hq that appears in the 
standard FRW model and measured in kms 'Mpc Setting 
flo = fl(fo) = 1 and Hq = H{to), we are able to write So and D as 
functions of s and Ho. Since the effective gravitational coupling 



is 



- TF , it turns out that in order to recover an attractive 



gravity we get s e (-2 , -1). Restricting furthermore the values 
of s to the range s e (-| , -1) the potential V(<p) is an inverse 
power-law, (p^^^P^"^^. In this framework, we obtain naturally an 
effective cosmological constant: 



2s + 3 



*^9) A.// = GeffP<p, 



(21) 
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where lim,^oo J^eff - Am + 0. It is then possible to associate an 
effective density parameter to the A term, via the usual relation 



A, 



eff 



(22) 



Eqs. (fT2l l and ( fT3l l are all that is needed to perform the t - z 
analysis. It is worth noting that, since the lookback time - red- 
shift (f - z) relation does not depend on the actual value of fo, 
it furnishes an independent cosmological test through age mea- 
surements, especially when it is applied to old objects at high 
redshifts. For varying w, as in the case of our model, the look- 
back time can be rewritten in a more general form by considering 



Jo 



Hod a 



aH(a) 
and then writing 



tdz) 



r Hoda 
Jo aH{a) ' 



where a = 1/(1 + z) and f// = ^ is the Hubble time. 



2.0.2. The case of quartic potentials 

As it is clear from Eq.(fT2ll and Eq.(fT3l) for a generic value of s 
both the scale factor a{t) and the scalar field (pit) have a power 
law dependence on time. It is also clear that there are some 
additional particular values of s, namely s - Q and 5 = -3 
which should be treated independently. Actually s - re- 
duces to the minimally coupled case (see |Demianski et al., 2005| l 
while s - -3) corresponds to the induced gravity with a quar- 
tic potentials: namely, F - -^(f^, and V{<p) - A(j)^ . This case 
is particularly relevant since it allows one to recover the self- 
interaction potential term of several finite temperature field the- 
ories. In fact, the quartic form of potential is required in or- 
der to implement the symmetry restoration in several Grand 
Unified Theories. Consequently, we limit our analysis to this 
model, showing that it naturally provides an accelerated ex- 
pansion of the universe and other interesting features of dark 
energy models. It worth to note that since the special deriva- 
tion of our model, we can focus the analysis on quantities 
concerning, mostly, the background evolution of the universe, 
without recurring to any observable related to the evolution of 
perturbations. The following analysis can be easily extended 
to more general classes of non-minimally coupled theories, 
where exact solutions can be achieved by Noether symmetries 
( [Capozziello et al, 1996| |Capozziello & de Ritis, 199 3 ), but, for 
the sake of simplicity, we restrict only to the above relevant case. 
The general solution is 



a{t) = aie— [(-1 + e°" ') +a2t + aiY , 



^4 



(-l+e"'i ') +02 z+ffi ' 



(23) 
(24) 



where ai = 4 VI, and a,- , 02 and 0,- are integration constants, 
related to the initial matter density, D, by the relation a\a2 - D, 
which implies that they cannot be null. The case A = has to be 
treated separately. 

It turns out that aj,, (pi and a,- have an immediate physical 
interpretation: a, is connected to the value of the scale factor at 
t - tQ - \, while Qf3 set the value fl(0). We can 

safely set a{Q) - 0, so that a-}, - 0. Actually this is not strictly 
correct, as our model does not extend up to the initial singularity. 



However, this position introduces a shift in the time scale which 
is small with respect to that of the radiation dominated era. This 
alters by a negligible amount the value of fo, while, moreover, it 
is left undetermined in our parametrisation. 

Since F(0) oc 0^, and Gf//(0) oc -j^, we note that a way 
of recovering attractive gravity is to chose 0,- as a pure imaginary 
number Without compromising the general nature of the prob- 
lem, we can set 0, = (, so that our field in Eq. ( |24] | becomes a 
purely imaginary field, giving rise to an apparent inconsistency 
with the choice in the Lagrangian (|2]l. Actually, the generic in- 
finitesimal generator of the Noether symmetry is 



X- a— +13— +a— +p— 

da d(f> da 50 ' 

where a and /3 are both functions of a and (p, and: 

da da . . d/3 d/3 . 
a = —a + —(f> ; p = —a + —(/>■ 
da d(f> da d(p 



(25) 



(26) 



Demanding the existence of Noether symmetry X.xL = 0, we get 
the following equations. 



a + 2a 



da jdBF' F' „ 
da da F F 



da dB\ , „ ^ da a^ dB ^ 
2a + a— +a^\F' +aF"/3 + 2F— + — = 
da d(f>l d(f) 6 da 



3a+l2F'(ip)—+2a^=Q 



r F' 
— = p{s) — 
V ' F 



(27) 



(28) 



(29) 



(30) 



It turns out that these equations are preserved if we adopt the 
more common signature (-, -H, -H, -h) in the metric tensor and flip 
the sign of the kinetic term. This imply that the generic infinites- 
imal generator of the Noether symmetry is preserved, and hence 
provides phantom solutions of the field equations. Therefore, as 
far as the quaestio of setting 0,- - i concerns, it turns out that the 
scalar field can be physically interpreted as a phantom solution 
of the field equations adopting the signature (-, +, -H, -h), even 
if it is mathematically represented as a pure imaginary field. It 
worth noting that such phantom solutions in scalar tensor gravity 
do not produce any violation in energy conditions, since p^, + 
p^ an d p^+3p^ are strict ly positive in the minimally coupled case 
(see |Ellis & Madsen, 1991| and references therein for a discus- 
sion). Furthermore, complex scalar fields (in particular purely 
imaginary ones) have been widely considered in classical and 
quantum cosmology with minimal and non-minimal couplings 
giving rise to interesting boundary conditions for inflationary be- 
haviors [Amendola et al., 1994| [Kamenshchik et al., 1995| 
Finally, assuming H(l) - 'Hq, we get 



a{t) 



where 



(3'Wo(c°i-l)-c?i(c''i+l))f 
yHa+a,-2 



1} 



2/3 



02(f) = - 



(3Wo(.°i-i)-»i(»°Ui)> I ' 

' 3'«„+ff,-2 ' 



3"Ho + ai-2 
2 + 2e"^{ai - 1) 



2/3 



(31) 



(32) 



(33) 
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We have to note that 6 = 



^eff 



39{q - a. 



and still ai 



Geff 2 

3'Ko, as we shall see in the next section, this expression is very 
small, as indicated by observations. This fact means that all 
physical processes implying the effective gravitational constant 
are not dramatically affected in the framework of our model. 

2. 1 . Some remarks 

Before starting the detailed description of our method to use the 
age measurements of a given cosmic clock to get cosmological 
constraints, it worth considering some caveats connected with 
our model. We concentrate in particular on the Newtonian limit 
and Post Parametrized Newtonian (PPN) parameters constraints 
of this theory, and discuss some advantages related to the use of 
the t(z) relation, with respect to the magnitude - redshift one, to 
constrain the cosmological parameters. 

2.1.1. Newtonian limit and Post Parametrized Newtonian 
(PPN) beliaviour 

Recently the cosmological relevance of extended gravity theo- 
ries, as scalar tensor or high order theories, has been widely 
explored. However, in the weaklimit approximation, all these 
classes of theories should be expected to reproduce the Einstein 
general relativity which, in any case, is experimentally tested 
only in this limit. This fact is matter of debate since several 
relativistic theories do not reproduce Einstein results at the 
Newtonian approximation but, in some sense, generalize them, 
giving rise, for example, to Yukawa like corrections to the 
Newtonian potential which could have interesting physical con- 
sequences. In this section we want to discuss the Newtonian limit 
of our model, and study the Post Parametrized Newtonian (PPN) 
behaviour. In order to recover the Newtonian limit, the metric 
tensor has to be decomposed as 



(34) 



where r/f^y is the Minkoskwi metric and /i^v is a small correction 
to it. In the same way, we define the scalar field i/r as a perturba- 
tion, of the same order of the components of hf,y, of the original 
field (A, that is 



(p = (f () + ((/, 



(35) 



where is a constant of order unit. It is clear that for - 1 
and 1/^ = Einstein general relativity is recovered. To write in an 
appropriate form the Einstein tensor G^y, we define the auxiliary 
fields 



and 



(36) 



(37) 



Given these definitions it turns out that the weak field limit of 
the power-law potential gives (see |Demianski et al., 2007| ) 



1 - 16^(^) 



M 



2^is)il - 12^(5)) 

Vo(p(s)-4)(p(s)-l) 



1 - n^{s) 



l-2^(^) 



(38) 



- An 



hij - 6ij 



2+p{s) 



Vo(p(s) - 4)(p(s) - I) 



2ipis)-l) 
1 - 8^(^) 



1-2^(5) 



M 
r 



Mfl Vo(p(s)-4)(p(s)-l) 



+ 4-7t 



M Voipis)-4)ipis)-l) 



(39) 



2(p{s) - 1) 



l-2^(s) 



where only linear terms in Vo are given and we omitted the con- 
stant terms. 

In the case of the quartic potential case, instead we obtain: 



100 



[1^5(1- 


16 a 


M 




[2 ^(1- 


120 \ 


r 


[ ^ \ 



hii ^ 6ii 



1 - 8^) 



1 



2^1-12^) 
2 



M 



0, 



+0 



2(1-12^)^0 



M 

— + (Ao, 



(40) 
(41) 
(42) 



A and 



being = 2(f^^\j/Q a sort of cosmological term, ^ 
iffQ an arbitrary integration constant. As we see, the role of 
the self-interaction potential is essential to obtain corrections 
to the Newtonian potential, which are constant or quadratic 
as for other generalized theories of gravity. Moreover, in gen- 
eral, any relativistic theory of gravitation can yield corrections 
to the Newton potential (see for example Will, 1993) which in 
the post-Newtonian (PPN) formalism, could furnish tests for 
the same theory based on local experiments. A satisfactory de- 
scription of PPN limit for scalar tensor theories has been devel- 
oped in (Damour & Esposito-Farese, 1993 1. The starting point 



of such an analysis consist in redefine the non minimally coupled 
Lagrangian action in term of a minimally coupled scalar field 
model via a conformal transformation from the Jordan to the 
Einstein frame. In the Einstein frame deviations from General 
Relativity can be characterized through Solar System experi- 
ments ( Will, 1 993] ) and binary pulsar observations which give 
an experimental estimate of the PPN parameters, introduced by 
Eddington to better determine the eventual deviation from the 
standard prediction of General Relativity, expanding the local 
metric as the Schwarzschild one, to higher order terms. The gen- 
eralization of this quantities to scalar-tensor theories allows the 
PPN-parameters to be expressed in term of the non-minimal cou- 
pling function F{(p) : 



/™ - 1 = - 



- 1 = 



F(<f>) + 2[F'(cf,)r 

F(4>yF'(4>) dy''™ 
F(cf,) + 3[F'{<p)]^ d<p 



(43) 



(44) 



Results about PPN parameters are summarized in Tab[T] The ex- 
perimental results can be substantially resumed into the two lim- 
its (jSchimJetaUjO^ 

lro™-l|<2xIO-\ lySf^- l|<6xl0-^ (45) 

If we apply the formulae in the Eqs. (|43]) and ( l44b to the power 
law potential, we obtain : 



1 



4^(5) 
1 + 8^*) 



(46) 
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Mercury Perih. Shift 


127™^ 11 < 3x10-^ 


Lunar Laser Rang. 


<™-rr^-3 = -(0.7±l)xl0-* 


Very Long Bas. Int. 


|7™-1| = 4x10-4 


Cassini spacecraft 


- 1 = (2.1±2.3)xl0-' 



Table 1. A schematic resume of recent constraints on the PPN- 
parameters from Solar System experiments. 



- 1 = 



PPN 



F(cf>yF'{cf>) dr 

F(0) + 3[F'(0)]2 d(f> 



= 0. 



(47) 



The above definitions imply that the PPN-parameters in gen- 
eral depend on the non-minimal coupling function F(<p) and 
its derivatives. However in our model y^™ depends only on 
s while y6^™ - Lit turned out that the limits for in 
the Eq. (l45T l is naturally verified, for each value of s, while 
the constrain on ly^™ ~ 1 | satisfied for s e (-1.5, -L4) 
(see IDemianski et al., 20 07V The case of the quartic potential 
is quite different: we easily realize that the constrain on 
is not satisfied, while it is automatically satisfied the constrain 
on 13^™, as well the further limit on the two PPN-parameters 
y^™ and /3^™, which can be out lined by means of the ratio 
( Damour & Esposito-Farese, 1993[ l: 



1 



< 1.1 . 



(48) 



However the whole procedure of testing extended theories of 
gravity based on local experiments (as PPN parameter con- 
straints) subtends two crucial questions of theoretical nature: 
the question of the so called inhomogeneous gravity, and the 
question of the conformal transformations. The former mainly 
consist in matching local scales with the cosmological back- 
ground: there is actually no reason a priori why local ex- 
periments should match behaviours occurring at cosmological 
scales. A non-minimally coupling function F{^), for instance, 
can alter the Hubble length at equivalence epoch, which is a scale 
imprinted on the power spectrum. CMBR, large-scale structure 
and, in general, cosmological experiments could provide com- 
plementary constraints on extended theories. In this sense also 
the coupUng parameter ^ oc ^ may be larger than locally is 
( [Clifton et al., 2004| [Acquaviva et al., 2005) . It is actually very 
interesting to compare just Umits on the Brans-Dicke parameter 
ojbd coming both from solar system experiments, (Dbd > 40000 
( |Will,1993| ), and from cuiTent cosmological observations, in- 
cluding cosmic microwave anisotropy data and the galaxy power 
spectrum data, cobd > 120) ( jAcquaviva et al., 2005 1. As fur- 
ther consideration, it is worth stressing that our model F{(p) oc 
(0 + (po)^ has been derived by requiring the existence of Noether 
symmetries for the pointlike Lagrangian defined in the FRW 
minisuperspace variables {a, a, 0,0): the extrapolation of such 
a solution at local scales is, in principle, not valid since the 
symmetries characterizing the cosmological model are not work- 
ing. However, these considerations do not imply that local ex- 
periments lose their validity in probing alternative theories of 
gravity [Will, 1 993[ but point out the urgency to understand how 
matching observational results at different scales in a coherent 
and self-consistent theory not available yet. Also the question re- 
lated to the conformal transformations is of great interest: these 



' In cosmology, we assume the Cosmological Principle and then a 
dynamical behavior averaged on large scales. This argument cannot be 
extrapolated to local scales, in particular to Solar System scales, since 
anisotropies and inhomogeneities cannot be neglected 



ones are often used to reduce non minimally coupled scalar field 
models to the cases of minimally coupled field models gain- 
ing wide mathematical simplifications. The Jordan frame, in 
which the scalar field is nonminimally coupled to the Ricci cur- 
vature, is then mapped into the Einstein frame in which the 
transformed scalar field is minimally coupled to the Ricci cur- 
vature, but where appears an interaction term between standard 
matter and the scalar field. It turns out the two frames are not 
physically equivalent, and some care have to be taken in apply- 
ing such techniques (see for instance (Faraoni, 2000] ) for a criti- 
cal discussion about this point). Actually even if a scalar tensor 
theory of gravity in the Jordan frame can be mimicked by an 
interaction between dark matter and dark energy in the Einstein 
frame, the influence of gravity is quite different in each of them. 
These considerations, however, do not imply that local experi- 
ments lose their validity to probe scalar tensor theories of grav- 
ity, but should simply highlight the necessity to compare local 
and cosmological results, in order to understand how to match 
different scales. 



2.1 .2. Comparing between the t(z) and the (m - M)(z) 
relations 

It is well known that in scalar tensor theories of gravity, as well in 
general relativity, the expansion history of the universe is driven 
by the function //(z); this implies that observational quantities, 
like the luminosity distance, the angular diameter distance, and 
the lookback time, are all function of H{z) (in particular H{a) ac- 
tually appears in the kernel of ones integral relations). It turns out 
that the most appropriate mathematical tool to study the sensi- 
tivity to the cosmological model of such observables consists in 
performing the functional derivative with respect to the cosmo- 
logical parameters (see [Saini et al, 2003[ for a discussion about 
this point in relation to distance measurements). In this section 
however we face the question from an empirical point of view: 
we actually limit to show that the lookback time is much more 
sensitive to the cosmological model than other observables, like 
luminosity distances, and the modulus of distance. This circum- 
stance encourages to use, together with other more standard 
techniques (as for instance the Hubble diagram from Snela ob- 
servations), the age of cosmic clocks to test cosmological scenar- 
ios, which could justify and explain the accelerated expansion of 
the universe. 

In any case, before describing the lookback method in de- 
tails, let us briefly sketch some distance-based methods in order 
to compare the two approaches. It is well known that the use of 
astrophysical standard candles provides a fundamental mean of 
measuring the cosmological parameters. Type la supernovae are 
the best candidates for this aim since they can be accurately cali- 
brated and can be detected up to enough high red-shift. This fact 
allows to discriminate among cosmological models. To this aim, 
one can fit a given model to the observed magnitude - redshift 
relation, conveniently expressed as : 



I2(z) = 5 log —dL(z) + 25 

no 



(49) 



being jj the distance modulus and di(z) the dimensionless lumi- 
nosity distance. Thus d^iz) is simply given as : 

ddz) = (1 + z) r dz'[^„,(l + z'f + Qa]-'^' . (50) 
Jo 

where Qa = 1 - Q„, in the case of ACDM model, but, in general, 
Qa can represent any dark energy density parameter. 
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The distance modulus can be obtained from observations of 
SNela. The apparent magnitude m is indeed measured, while 
the absolute magnitude M may be deduced from template fitting 
or using the Multi- Color Lightcurve Shape (MLCS) method. 
The distance modulus is then simply ji - M - m. Finally, the 
redshift z of the supernova can be determined accurately from 
the host galaxy spectrum or (with a larger uncertainty) from the 
supernova spectrum. 

Roughly speaking, a given model can be fully characterized 
by two parameters : the today Hubble constant Hq and the matter 
density Q,„. Their best fit values can be obtained by minimizing 
the defined as : 



[//f°^fci//o,a„)- /if^]^ 

Alo,' inz,i 



(51) 



where the sum is over the data points. In Eg. (Bill. cr^„ is the es- 
timated error of the distance modulus and cr„,-. is the dispersion 
in the distance modulus due to the uncertainty cr, on the SN red- 
shift. We have : 



1 ddL 
In 10 \iiZ"5F 



(52) 



where we can assume cr^ - 200 km s ' adding in quadrature 
2500 km s"' for those SNe whose redshift is determined from 
broad features. Note that cr,„. depends on the cosmological pa- 
rameters so that an iterative procedure to find the best fit values 
can be assumed. 

For example, the High-z team and the Supernova 
Cosmology Project have detected a quite large sample of high 
redshift (z ^ 0. 18 - 0.83) SNe la, while the Calan - Tololo survey 
has investigated the nearby sources. Using the data in Perlmutter 
et al. [Perlmutter et al., 1997| and Riess et al. [Riess, 20^01 com- 
bined samples of SNe can be compiled giving confidence regions 
in the (Qm, ^^o) plane. 

Besides the above method, the Hubble constant Hq and the 
parameter Qa can be constrained also by the angular diameter 
distance Da as measured using the Sunyaev-Zeldovich effect 
(SZE) and the thermal bremsstrahlung (X-ray brightness data) 
for galaxy clusters. Distances measurements using SZE and X- 
ray emission from the intracluster medium are based on the fact 
that these processes depend on different combinations of some 
parameter of the clusters. The SZE is a result of the inverse 
Compton scattering of the CMB photons of hot electrons of the 
intracluster gas. The photon number is preserved, but photons 
gain energy and thus a decrement of the temperature is generated 
in the Rayleigh- Jeans part of the black-body spectrum while an 
increment rises up in the Wien region. The analysis can be lim- 
ited to the so called thermal or static SZE, which is present in 
all the clusters, neglecting the kinematic effect, which is present 
in those clusters with a nonzero peculiar velocity with respect to 
the Hubble flow along the line of sight. Typically, the thermal 
SZE is an order of magnitude larger than the kinematic one. The 
shift of temperature is: 



To 



= 3' 



jccoth 



(53) 



where x ■ 



hv 



is a dimensionless variable, T is the radiation 



temperature, and y is the so called Compton parameter, defined 
as the optical depth t - ctt J ngdl times the energy gain per 
scattering: 



CKbT, 

J m^c^ 



o-rdl. 



(54) 



In the Eq. ( |54] i, Tg is the temperature of the electrons in the intr- 
acluster gas, mg is the electron mass, «g is the numerical density 
of the electrons, and ctj is the cross section of Thompson elec- 
tron scattering. The condition » T can be adopted {7^ is the 
order of 10^ K and T, which is the CBR temperature is ^ 2.7 K). 
Considering the low frequency regime of the Rayleigh-Jeans ap- 
proximation one obtains 



To 



(55) 



The next step to quantify the SZE decrement is to specify the 
models for the intracluster electron density and temperature 
distribution. The most commonly used model is the so called 
isothermal /3 model. One has 



Heir) = Heir) 



1+1^ 

re. 



2\ 



(56) 
(57) 



being rte^ and T^g, respectively the central electron number den- 
sity and temperature of the intracluster electron gas, r^ and /3 are 
fitting parameters connected with the model. Then we have 



AT _ 
being 



2KBO-TTe„ tie 



/-.CO / / \ 2\ 2 



(58) 



(59) 



The integral in Eq. ( |59] l is overestimated since clusters have a 
finite radius. 

A simple geometrical argument converts the integral in 
Eq. ( |59l ) in angular form, by introducing the angular diameter 
distance, so that 



I = 0, 



2\ 



1/2-3^/2 



r(¥-| ) 
r(f) 



■ roR- 



(60) 



In terms of the dimensionless angular diameter distances, c/a 



(such that Da - — dA), one gets 

^0 



M{9) 



2 a-jKBTecrie 
Ho mpC 



r(f) 



02 



2\ 



i(I-3/S) 



dA,m) 



and, consequently, for the central temperature decrement, we ob- 
tain 



AT{6 = 0) 
7^ 



2 O-jKbT eerie 

Ho meC 



r(f)- 



Ho 



dA- 



(62) 



The factor — dA in Eq. ( l62l ) carries the dependence on the ther- 
^0 

mal SZE on both the cosmological models (through Ho and the 
Dyer-Roeder distance dA) and the redshift (through t/^)- From 
Eq. ( l62l ). we also note that the central electron number density 
is proportional to the inverse of the angular diameter distance, 
when calculated through SZE measurements. This circumstance 
allows to determine the distance of cluster, and then the Hubble 
constant, by the measurements of its thermal SZE and its X-ray 
emission. 
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This possibility is based on the different power laws, accord- 
ing to which the decrement of the temperature in the SZE, 
, and the X-ray emissivity, Sx, scale with respect to the electron 
density. In fact, as pointed out, the electron density, when calcu- 
lated from SZE data, scales as ( n^,f^ oc c/^'), while the same 



one scales as d^^ (n^g 



x-ray 



A ^ 'eO "A 

^ ) when calculated from X-ray 
data. Actually, for the X-ray surface brightness, S x, assuming 
for the temperature distribution of Tp = T^o, one gets the follow- 
ing formula: 



fx 2 1 , 



c 

•—dAhx, 
no 



(63) 



being 



Jo \neol 



dl 



the X-ray structure integral, and Sx the spectral emissivity of the 
gas (which, for > 3x10^, can be approximated by a typical 
value: sx - sy/T\, , with s ^ 3.0xlO"^^n^ erg cm^^ i"' ■^^"')- 
The angular diameter distance can be deduced by eliminating the 
electron density from Eqs. (l62T i and (1631 ), yielding: 



^ 47T(l+z)\ JkBcrT\ 
Sx e 



3/2, 



. — dAX- 

//() B(3/3- 



1 \ 

2' 2 



1^ 

2' 2 



(64) 



T(a)r(b) 

where B(a, b) - is the Beta function. 

T{a + b) 

It turns out that 



c , (AT-q)^ 1 

= dA 

Ho SxoTl « 



(65) 



eO 



where all these quantities are evaluated along the line of sight 
towards the center of the cluster (subscript 0), and ft- is referred 
to a characteristic scale of the cluster along the line of sight. It 
is evident that the specific meaning of this scale depends on the 
density model adopted for clusters. In general, the so called /3 
model is used. 

Eqs. ( l64b allows to compute the Hubble constant Hq, once 
the redshift of the cluster is known and the other cosmological 
parameters are, in same way, constrained. Since the dimension- 
less Dyer-Roeder distance, d^, depends on Qa, Q,„, comparing 
the estimated values with the theoretical formulas for D^, it is 
possible to obtain information about 0,„, Oa, and Hq. Modelling 
the intracluster gas as a spherical isothermal yS-model allows 
to obtain constraints on the Hubble constant Ho in a standard 
ACDM model. In general, the results are in good agreement with 
those derived from SNe la data. 

Apart from the advantage to provide an further instrument of 
investigation, the lookback time method uses the stronger sen- 
sitivity to the cosmological model, characterizing the t{z) rela- 
tion, as shown in Figs. ( I3I6I I. Moreover, as we will discuss in the 
next sections such a method reveals its full validity when applied 
to old objects at very high z: actually it turns out that this kind 
of analysis is very strict and could remove, or at least reduce, 
the degeneracy which we observe at lower redshifts, also, for 
instance, considering the Hubble diagram observations, where 
different cosmological models allow to fit with the same statis- 
tical significance several observational data. In a certain sense, 
we could argue that the lookback time, even if exhibits a wide 
mathematical homogeneity with the distance observables, does 



not contain the same information, but, rather, presents some in- 
teresting peculiar properties, useful to investigate also alterna- 
tive gravity theories. However it is important to remark that the 
comparison with observational data in scalar tensor theories of 
gravity is more complex than in general relativity, since the ac- 
tion of gravity is different. For instance, the use of type la super- 
novae to constraint the cosmological parameters (and hence the 
claim that our universe is accelerating) mostly lies in the fact that 
we believe that they are standard candles so that we can recon- 
struct the luminosity distance vs redshift relation and compare it 
with its theoretical value. In a scalar-tensor theory, we have to 
address both questions, i.e: the determination of the luminosity 
distance vs redshift relation and the property of standard candle 
since two supernovae at different redshift probe different grav- 
itational coupling constant and could not be standard candles 
anymore. Actually it can be shown that in this case it is possible 
to generalize the theoretical expression of the distance modulus, 
taking into account the effect of the variation of Gejf through 
a correction term (Gaztanaga et al., 2001 1. Of course, the varia- 



tion of Geff with time (and then redshift) could challenges the 
reliability of age measurements of some cosmic clocks to test 
cosmological models, unless to construct some reasonable theo- 
retical model, which quantifies and corrects this effect. Actually 
the estimation of cluster ages is based on stellar population syn- 
thesis models which rely on stellar evolution models formulated 
in a Newtonian framework. In order to be consistent with scalar 
tensor gravity, which assume an evolving gravitational coupling 
Geff, one should investigate to which extent this Ge/f variation 
affects the results of the given stellar population synthesis model. 
However, it turns out that in the context of our qualitative anal- 
ysis, the effective gravitational coupling varies in a range of no 
more than 6% and then the effect of such a variation can be in- 
cluded in a bias factor df, which, resulting to be df - 4.5 ± 0.5, 
as we will see in the next section, affects the age measurements 
of the considered cluster sample much more than the variation 
ofG,//. 



3. The dataset at low and intermediate redshifts 

In order to discuss age constraints for the above background 
models, we first use the dataset compiled by Capozziello et 



al. ^Capozziello et al., 2004 1, and given in Table 1, which con- 
sists of age estimates of galaxy clusters at six redshifts dis- 
tributed in the inte rval 0.10 < z < 1.27 (see Sec. IV of ( 
[Capozziello et al.,~2 004) for more details on this age sample). 
To extend such a dataset to higher redshifts, we join the GDDS 
sample presented by McCarthy et al. (2004), consisting of 20 old 
passive galaxies, distributed over the redshift interval 1.308 < 
z < 2.147, and shown in Fig. ([T]). In order to build up our total 
loockback time sample, we first select from GDDS observations 
what we will consider as the most appropriate data to our cos- 
mological analysis. Following (Dantas et al |Dantas et al., 2007| 
we adopt the criterion that given two objects at (approximately) 
the same z, the oldest one is always selected, ending up with a 
sample of 8 data points, as in Fig.©. 

Through the Eqs. dA.ll ) and ( IA.2I ) in Appendix A, which in 
our case can only be evaluated numerically, we perform a 
analysis. For the power law potential we obtain, in our units, 
= 0.95, 'Hq = l.OOtZv ' = -l-39!o:4, = 14.04 ± .08 
Gyr, and df = 3.6 + 0.7 Gyr. Such best fit values correspond to 
Q^^^^ = 0.76+g gg, according to the Eq. ( |22] |. 

In the case of quartic potential we obtain - 1.01, 'Hq = 
13.04 ± .05 Gyr, and df = 4.2 + 0.7 



1.00!»:;»,a3 



o c+0.5 , 
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Color age 


Scatter age 


z 


N 


Age (Gyr) 


z 


N 


Age (Gyr) 




0.60 


1 


4.53 


0.10 


55 


10.65 




0.70 


3 


3.93 


0.25 


103 


8.89 




0.80 


2 


3.41 


1.27 


1 


1.60 





1.2 1.4 1.6 l.i 



2.2 2.4 2.6 



Table 2. Main properties of the cluster sample compiled by 
Capozziello et al. (Capozziello et al., 2004) used for the analy- 
sis. The data in the left part of the Table refers to clusters whose 
age has been estimated from the color of the reddest galaxies 
(color age), while that of clusters in the right part has been ob- 
tained by the color scatter (scatter age). For each data point, we 
give the redshift z, the number of clusters used and the age 
estimate 



Fig. 1. Original data from GDDS. This sample corresponds to 20 old 
passive galaxies distributed over the redshift interval 1 .308 < z < 2. 147, 
as given by McCarthy et al. (2004). 




0.5 1 1.5 2 2.5 3 



Fig, 4, The observational data of the whole dataset fitted to our quartic 
potential model, with '74 = l.OO+SJ^^, 0-3 = 2.5+j;^, to = 13.04± .05 Gyr, 
and df = 4.2 ± 0.7 Gyr. 



Fig. 2. The 8 high-z measurements selected after the criterion discussed 
in the text. 



Gyr. In the Figs. (O and |4]i, the observational data are plotted vs 
our best fit cosmological model. 

Remark: The range of values for "Ko does not correspond to 
a variation in the physical value of //q, which is a prior for the 
model. It reflects instead the scatter in the universe age. 




Fig. 3. The observational data of the whole dataset fitted to our model. 



with 'Ho = 1.00: 
3.6 ± 0.7 Gyr. 



+0.01 

0.07 ' 



-1.39+ 



14.04 ± .08 Gyr, and df 



4. Extending the analysis at high redshifts 

Previous discussion shows that the scalar-tensor quintessence 
model which we are studying can be successfully constrained 
by cosmic clocks (clusters of galaxies) at low (z ~ H- 0.5) 
and intermediate (z ~ 1.0 1.5) redshifts. In this section we 
investigate its viability vs the age estimates of some high red- 
shift objects, with a minimal stellar age of 1.8 Gyr, 3.5 Gyr and 
4.0 Gyr, respectively. It is actually well known that the evo- 
lution of the universe age with redshift (-^) differs from a 
scenario to another; this means that models in which the uni- 
verse is old enough to explain the total expansion age at z = 
may not be compatible with age estimates of high redshifts 
objects. This reinforces the idea that dating of objects consti- 
tutes a powerful methods to constrain the age of the universe 
at different stages of its evolution (jKrauss & Chaboyer, 2003| 
IDunlop et al. 1996| [Ferreras & Silk, 2003| ), and the first epoch 
of the quasar formation can be a useful tool for discriminating 
among different scenarios of dark energy pimenez et al., 2003] 
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Fig. 5. We compare the sensitivity to the values of the parameters in 
our quartic potential model in the lookback time relation and in the 
modulus of distance. Actually we plot the relative variation in t(z) (up- 
per diagram) and in - M with respect to a variation of ai from 3. to 3.5 
(the other parameters are fixed). It turns out that the t(z) relation is much 
more sensible. 



0.2 




0.2 0.4 0.6 0.8 1 1.2 



z 



Fig. 6. We compare the sensitivity of the ^ relation to the values of the 
parameters in the power law potential cosmological model. Actually we 
plot the relative variation in ^ with respect to a variation of s from -1.4 

to -1.3 (red line), and with respect to a variation of Hq from 1 to .9 (blue 
line), the other parameters being fixed. 

[Alcanitz et al., 200 ?). The existence of some recently reported 
old high-redshift objects if of relevance here, namely LBDS 
53W091, a 3.5 Gyr-oldradio galaxy atz = 1.55, LBDS 53 W069, 
a 4 Gyr-old radio-galaxy at z = 1 .43, and APM 08279-1-5255, an 
old quasar at z = 3.91, whose age was firstly estimated to lie 
in the range 2-^3 Gyr (Hasinge r et al., 2002|, a nd then updated 
to the range 1.8 H- 2.1 Gyr ( |Friaga et al., 2005 j l. It is clear that 
these objects can be used to impose more strict constraints on 
our model. 

We take advantage for the fact that we have exact solutions, 
so that the redshift-time relation can be inverted. It is easily de- 
rived from Eq.([3TTl. 

Taking for granted that the universe age at any redshift must 
be greater than or at most equal to the age of the oldest object 
contained in it, we introduces the the ratio 

■??=-> 1, 

as in ( [Alcanitz et al., 2003] l, with t, derived from Eq. (ISTT i. and 
tg the measured age of the objects. For each extragalactic object. 



Fig. 7. The dashed line indicates T(z) = HQt(z) according to the best 
fit values of Sec. IV. for the power law potential model. The shadowed 
region covers the allowed ranges coming from lookback time tests. The 
vertical lines give the lowest ages Tgof the considered objects for the 
allowed range of h = (0.64 H- 0.80). The figures reported in the text 
correspond to /z = 0.64, i.e. the lowest edges of the lines. For the 
first two objects, there is full agreement with the requirement 
T > Tg. The third age is marginally compatible, likely due to 
systematics in evaluating the age-metalicity relation. 

this inequality defines a dimensionless parameter Tg = Hotg. 
In particular, for LBDS 53W091 radio galaxy discovered by 
Dunlop et al. ( [Dunlop et al. 1996 ), the lower limit for the age 
yields tg{z = 1.55) = 3.5i/o Gyr which takes values in the inter- 
val 0.21 < Tg < 0.28. It therefore follows that Tg > 0.21. Similar 
considerations may also be applied to the 4.0 Gyr old galaxy at 
z = 1.43, for which Tg > 0.24, and to the APM 08279-H5255, 
which corresponds to Tg > 0.131. 

Only models having an expanding age parameter larger than 
Tg at the above values of redshift will be compatible with the 
existence of such objects. 

In Figs. (|7]i, and dHJ we show the diagram of the dimen- 
sionless age parameter T{z) - Hotiz) as a function of the 
redshift for two extreme best fit values of the parameters (s 
and ^Hq in the case of power law models, a-^ and "Ho for the 
quartic potential model), coming from the lookback time test. 
We observe that the age test turns out to be quite critical, 
even though such values of the parameters are fully compat- 
ible and equivalent to the analysis of other datasets (see 
IDemianski et al., 2006|[Demianski et al., 2007| l. 

In fact, it turns out that the value Tg > 0.131 
( [Friaga et al., 2005"} at z = 3.91 is quite selective in both cases. 

A remark is necessary at this point: the range of ages for 
quasar APM 08279+5255 is 1.8 and 2 Gyr. It is required for 
models from 10" M© to 10'^ M© in order to reach the metalicity 
Fe/0=2.5x solar unit (for details see |Friaga et al., 2005} . 

5. Discussion and Conclusions 

Cosmological models can be constrained not only using distance 
indicators but also cosmic clocks, once efficient methods are de- 
veloped to estimate the age of distant objects. Among these, the 
relations between lookback time and redshift z are particularly 
useful to discriminate among the huge class of dark energy mod- 
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Fig. 8. The same than in Fig.([8ll for the quartic potentail model 

els which have been recently developed to explain the observed 
present day accelerated behaviour of cosmic flow. In a near fu- 
ture, beside distance measurements, time measurements could 
greatly contribute to achieve a final cosmological model con- 
straining the main cosmological parameters. 

In this context, it is of fundamental importance to obtain 
valid cosmic clocks at low, intermediate and high redshifts in 
order to fit, in principle, a model at every epoch. In this case, 
degeneracies are removed and the reliability of the model can be 
proved. 

In this paper, we have tested a non minimally coupled 
scalar-tensor quintessence model, characterised by quartic self- 
interacting potential, using, as cosmic clocks, cluster of galaxies 
at low and intermediate redshifts, two old radio galaxies at high 
redshift, and a very far quasar The results are comfortable since 
the model seems to work in all the regimes considered. However, 
to be completely reliable, the dataset should be further enlarged 
and the method considered also for other time indicators. 

The main feature of this approach is the fact that cosmic 
clocks are completely independent of each other, so that, in prin- 
ciple, it is possible to avoid biases due to primary, secondary and 
so on indicators, as in cosmic ladder method. In that case, every 
rung of the ladder is affected by the errors of former ones and it 
affects the successive ones. By using cosmic clocks, this short- 
coming can be, in principle, avoided, since indicators are, by 
definition, independent. In fact, we have used different methods 
to test the model at low and high redshift with different indi- 
cators, which seem to confirm independently the proposed dark 
energy model. 

Another comment is due at this point. Having normalized 
the Hubble parameter at present epoch, assuming 'Ko = ^^(1), 
we do not need, in principle, priors on such a parameter, since 
we are using an exact solution. In other words, we can check the 
validity of the model by selecting reliable cosmic clocks only. 
Another advantage of such a choice is that one has to handle 
only small numbers of parameters in numerical computations. 
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Appendix A: The lookback time method 

In order to use the age measurements of a given cosmic clock 
to get cosmological constraints, let us consider an object / at 
redshift z and let f,(z) be its age defined as the difference between 
the age of the universe at the formation redshift zf, and at z: 



tiiz) = tdzp) - tdz). 

If one is able to to estimate the ages {?,) for / - 1,2,.. 
objects, we can estimate the lookback time tf'^izd as 



tTizd = tdZF) - tiiz) 

ti(z) - df , 



(A.l) 
,NofN 



(A.2) 



where f° ' is the age of the universe (which in our units is set to 
1), while the bias (delay factor) can be defined as : 



df^t\ 



.obs 



tliZf) 



(A3) 



The delay factor is introduced to take into account our ignorance 
of the formation redshift zf of the object. Actually, to estimate 
Zf, one should use Eq. dA.ll i, assuming a background cosmolog- 
ical model. Since our aim is to constrain the background cos- 
mological model, it is clear that we cannot infer zf from the 
measured age, so that this quantity is a priori undetermined. 
Moreover we rely that it can take into account also the effect of 
the Geff variation on the age estimations, because the expected 
magnitude of such an effect. In principle, df should be different 
for each object in the sample unless there is a theoretical reason 
to assume the same redshift at the formation of all the objects. 
However, we can realistically assume that df is the same for all 
the homologous objects of a given dataset (in the range of the 
errors), and we consider df rather than zf as the unknown pa- 
rameter to be determined from the data. 
We may then define a merit function : 



with h°^-^ the estimated value of h and cr^ its uncertainty. We 
use the HST Key project results (Freedman et al., 2001[ ) setting 
(h, cTfi) - (0.72, 0.08). Note that this estimate is independent of 
the cosmological model since it has been obtained from local 
distance ladder methods. 

The best fit model parameters p may be obtained by maxi- 
mizing -C(p) which is equivalent to minimize the defined in 
Eq. dA.Sb . It is worth stressing that such a function should not be 
considered as a statistical in the sense that it is not forced to be 
of order 1 for the best fit model to be considered as a successful 
fit. Actually, such an interpretation is not possible since the er- 
rors on the measured quantities (both f,- and to) are not Gaussian 
distributed. Moreover, there are uncontrolled systematic uncer- 
tainties that may also dominate the error budget. Moreover, there 
are uncontrolled systematic uncertainties that may also dominate 
the error budget. Nevertheless, a qualitative comparison of dif- 
ferent models may be obtained by comparing the values of this 
pseudo x^, even if this should not be considered as definitive 
evidence against a given model. 

Given that we have more than one parameter, we obtain the 
best fit value of each single parameter as the value which max- 
imizes the marginalized likelihood for that parameter, defined 
as : 

Xp, oc J dpi . . . J dpi-i J dpi+i . . . J dp„ £(p) . (A.6) 

After having normalized the marginalized likelihood to 1 at max- 
imum, we compute the 1 cr and 2cr confidence limits (CL) on that 
parameter by solving Xp, - exp(-0.5) and X.p. = exp(-2) re- 
spectively. 



1 



Np + 1 



i-theor 



(P)- 



^obs 



z 



t'r'(zi,p)-tr(zi) 



■2 + (r2 



where A^^ is the number of parameters of our model, cr,, cr, are 
the uncertainties on fjj''' and f^Hzi). Here the superscript theor 
denotes the predicted values of a given quantity. 

In principle, such a method can work efficiently to discrim- 
inate between the various cosmological models. However, the 
main difficulty is due to the lack of available data which leads 
to large uncertainties on the estimated parameters. In order to 
partially alleviate this problem, we can add further constraints 
on the model by using priors; for example choosing a Gaussian 
prior on the Hubble constant allows us to redefining the likeli- 
hood function as 



-£(p) Z„(^) exp 



1 Ih-h"'" 
2 1 o-h 



(xexp[-;^2(p)/2], (A.4) 



where we have absorbed df lo the set of parameters p and have 
defined : 

uobs \ ^ 



X 



--xl + 



'h-h" 

\ o-h 



(A.5) 



